Universality of short-range nucleon-nucleon correlations 
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Short-range correlations between nucleon pairs in different spin-isospin channels are investigated for light 
nuclei using the Argonne v8' interaction. At distances below 1 fm a universal behavior is found for the deuteron, 
^H, ^He and for ground and first excited states in ^He. This behavior in coordinate space is reflected by a 
universal behavior for the high-momentum components in momentum space. The universality indicates that a 
pairwise renormalization is possible in order to obtain a universal eff'ective two-body interaction that does not 
scatter to high momentum states. The exact two-body densities are compared with those obtained using the 
unitary correlation operator method with simple trial wave functions. The eff'ect of three-body correlations due 
to the tensor force on the two-body densities is discussed. 
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I. INTRODUCTION 

Realistic nucleon-nucleon interactions, which reproduce 
the scattering phase shifts, imply usually strong repulsive and 
tensor forces at short distances. These induce short-range cor- 
relations in the nuclear many-body system which complicate 
the theoretical description so much that exact solutions of the 
many-body Schrodinger equation become unfeasible for sys- 
tems with more than about twelve nucleons. Therefore theo- 
retical methods have to be devised in order to tackle this prob- 
lem. 

At short distances, where the scattering nucleons overlap 
strongly, there is no unique way to parametrize the complex 
many-body quantum chromodynamics problem in terms of 
just the distance, the relative momentum, and the spins of the 
nucleons. In all models for the nucleon-nucleon interaction 
the short-range behavior is governed by form factors of vari- 
ous types without rigorous derivations. Therefore experimen- 
tal data for elastic scattering, which provide phase shifts for 
energies up to the pion threshold cannot sufficiently constrain 
the nucleon-nucleon potential at small distances. In conse- 
quence, different phase shift equivalent interactions, (e.g., [1- 
4]), show a quite different high-momentum or short-range be- 
havior. 

Another uncertainty arises when going from the two- 
nucleon scattering states to bound many-body states of nuclei. 
In the scattering situation the two nucleons are in an energy 
eigenstate with a well-defined relation between momentum or 
kinetic energy, potential energy, and total energy, usually la- 
beled with on- shell. In the many-body case a nucleon pair, 
which interacts with the surrounding other nucleons, neither 
has sharp energy nor is their relative momentum related to 
their energy. In this situation the so-called off- shell behavior 
(i.e., local versus momentum-dependent parts) of the nuclear 
interaction is important but also not fully constrained by scat- 
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tering data. This ambiguity in the off- shell behavior of two- 
body interactions is also related to the three-body interactions 
that should accompany the different two-body interactions [5]. 

Information about the short-range behavior of the nuclear 
interaction is contained in the one- and two-body momentum 
distributions of nucleons in finite nuclei. However, it is dif- 
ficult to relate measured cross sections with momentum dis- 
tributions. Reactions and their kinematic conditions have to 
be chosen properly, for example, in such a way that the cor- 
related nucleon is removed instantaneously and final- state in- 
teractions are minimized [6, 7]. In recent years we have seen 
a renewed interest in studying short-range correlations using 
(e, e'pp) and (^, e' pn) [8, 9] (p,pp) and (p,ppn) [10] experi- 
ments. One important result from these studies is that the mo- 
mentum distributions above the Fermi momentum are domi- 
nated by tensor correlations [11]. There are also attempts to 
explore the effect of tensor correlations in nuclei by pick-up 
in (p, d) reactions [12]. 

Over the years short-range correlations have been studied 
using approaches such as the coupled-cluster method [13], 
correlated basis functions [14, 15], Green's functions [16], 
variational methods [17] or within a cluster expansion ap- 
proach [18, 19]. For a review see Ref. [20]. These meth- 
ods were essentially applied only to doubly-close shell nuclei 
like ^^O and "^^Ca. For lighter systems pioneering studies have 
been performed in the Green's function Monte Carlo approach 
[21-23] employing two- and three-body interactions. 

In this paper we do not investigate different realistic interac- 
tions but concentrate on the Argonne v8' (AV8') potential [24] 
where the short-range physics is described by a phenomeno- 
logical local potential. Extending on the results obtained in 
Ref. [25] the aim of this investigation is to solve the three- 
and four-nucleon system exactly and analyze the short-range 
correlations in the different spin-isospin channels. After ex- 
plaining briefly the many-body approach in Sec. II A we de- 
fine explicitly various one- and two-body densities that are 
used in Sec. Ill to illustrate that the short-range pair correla- 
tions are universal in the sense that they depend very little on 
the surrounding nucleons and the many-body state in general. 



This feature has been reaUzed some time ago by Forest et al 
[21]. Here we investigate the universaHty and discuss the im- 
pHcations to devise effective low-momentum interactions. 

The universal behavior gives hope that it is possible to de- 
rive effective interactions that are phase- shift equivalent and 
soft enough to permit many-body calculations with a Slater 
determinant basis without inducing large many-body effective 
potentials. The transformation of the Hamiltonian to an ef- 
fective one implies of course also the same transformation of 
any other observable. In order to keep the physics transparent 
it is highly desirable that these operators, which are usually 
one-body operators, do not take on large many-body contri- 
butions when they are transformed to effective ones. This will 
be the case if the induced correlations are of short range and 
the theoretical treatment takes account of this, as in the unitary 
correlation operator method (UCOM). Many observables like 
radii or multipole moments are not sensitive to short-range 
correlations, however, observables containing the spins, as in 
Gamow-Teller transitions are affected by the tensor correla- 
tions. 

In Sec. II we briefly recapitulate the many-body method 
adopted here, define the one- and two-body densities in coor- 
dinate and momentum space, and discuss the AV8' potential in 
the different spin-isospin channels. In Sec. Ill we display the 
correlations in coordinate and momentum space for the differ- 
ent spin-isospin-channels and investigate quantitatively their 
universality. In Sec. HIE we discuss how three-body corre- 
lations manifest themselves in the two-body densities and in 
Sec. Ill F we compare to two-body densities obtained with the 
UCOM approach. Summary and outlook are drawn in Sec. IV. 



II. MANY-BODY MODEL, INTERACTIONS, DENSITIES 
A. Correlated Gaussian basis approach 

We assume that an A-nucleon state can be expanded in 
terms of a combination of basis states, each of which is a prod- 
uct of space, spin and isospin parts. 



from the total center of mass of the system. The correlated 
Gaussian basis is conveniently expressed in terms of the rel- 
ative coordinates, e.g., the Jacobi set of the coordinates, x = 

(xu ^2, . . . ,^A-i): xi =ri- h. X2 = (h + h)/^ -h.--- The 
correlated Gaussian evaluated at the position x corresponding 
to the operator x takes the following form 
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where J/lm(^) = I^I^^lmC^) is a solid spherical harmon- 
ics. In Eq. (3) the matrix notation is used to simplify the ex- 
pression. That is, Ui is an A - 1 dimensional column vector 
and Ui denotes its transpose, uix = T/jll(ui)jXj. Similarly, Af 
is an (A - 1) X (A - 1) positive-definite, symmetric matrix, and 
xAiX is a short-hand notation for Yj%li(^i)jkXj • Xk. The basis 
is in fact correlated because all the coordinates are coupled 
through the off-diagonal elements of A/. The elements of A/ 
and Ui (and v/) are parameters to characterize the "shape" of 
the basis function. 

The matrix elements of the Hamiltonian can be analyti- 
cally obtained using the generating function technique. All 
the formulas needed are given in Refs. [25, 27, 28]. As 
seen in Eqs. (1), (2) and (3), each basis function contains 
both discrete and continuous parameters. The former includes 
Li/,L2/,L„5'i2,5'i23,...,5'/,ri2,ri23,...,(r/ = T) and the 
latter the elements of A/, ui, and v/. Though the Gaussians 
may not be ideal to cope with the repulsion, it is in fact possi- 
ble to obtain results as accurate as other sophisticated methods 
for a few-body system [27, 29]. One of the advantages of the 
present method is that the state T is expressed analytically so 
that physical quantities of interest can be readily evaluated. 
Since the Fourier transform of the correlated Gaussian basis 
is also expressed as correlated Gaussians in momentum vari- 
ables [27, 28], it is straightforward to calculate the matrix ele- 
ment of a quantity depending on the momentum operator. To 
have a compact basis size K saves time of computations. We 
use the stochastic variational method [26, 30] to choose the 
parameters and increase the basis dimension until good con- 
vergence is reached. 
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Here Ji is the antisymmetrizer and the square bracket [• • • ] 
stands for the angular momentum coupling. The spin and 
isospin parts are expanded using the basis of successive cou- 
pling, e.g.. 



B. One- and two-body density 

The antisymmetrized many-body state |^; JM) contains all 
the information about the nuclear system. For example, the 
one-body density in coordinate space is defined as 
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where the set of intermediate spins (Su.Sud,. . .) takes all 
possible values compatible with the total spin S i of the /th ba- 
sis. The isospin mixing is ignored in this paper, so that the 
total isospin Ti is kept fixed to T. The orbital part j^j^P^^^^ is 
expressed in terms of the explicitly correlated Gaussian ba- 
sis [26, 27] as explained below. 

We denote the position operator of particle / by r/. To 
simplify the notation, we define this position to be measured 
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where r/ is the position operator for the /th particle measured 
from the position of the total center of mass. Likewise the 
one-body momentum distribution is calculated as 
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where the momentum ki of particle / is defined in the total 
momentum frame of the nucleus. That means if a particle has 
momentum A: 1 all other particles have together a total momen- 
tum -ki. 

One should keep in mind that the possibility of finding a 
single nucleon with momentum ki does not imply that this 
nucleon has an energy that is related to ki, such as k\l{2mN). 
As all nucleons are interacting with each other one can not 
define an observable for the energy of one nucleon. 

In a mean field picture, where particles move independently 
in a common single-particle potential, each particle can be as- 
signed to a single-particle state that has a sharp energy, the 
single-particle energy. But this energy is also not uniquely 
related to a momentum, because the state has spread-out dis- 
tributions in momentum and in coordinate space which are 
related. 

Similar eff'ects occur in the interacting many-body case 
where only the total energy (i.e., the eigenvalue of the Hamil- 
tonian) is well defined. Rapid spatial variations in the many- 
body wave function show up as increased probabilities at large 
single-particle momenta. For example strongly repulsive two- 
body interactions induce areas where one will not find particle 
pairs because their interaction energy would be large and pos- 
itive. At the edges of these correlation holes the wave func- 
tion has to vary rapidly giving rise to large momenta and extra 
positive kinetic energy. But altogether it is more profitable 
to pay the kinetic energy and avoid the even larger poten- 
tial energy by staying out of the repulsive region. Thus the 
high-momentum tail of the momentum distribution reflects the 
short-range correlations. 

In Sec. Ill we discuss these phenomena by looking at dif- 
ferent exact many-body eigenstates. The one-body densities 
can be accessed in scattering experiments, the proton density 
preferably by electron scattering. 

A more direct way to see short-range correlations is given 
by the two-body density 



their center-of-mass position R = (ri -\-r2)/2. 
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The corresponding distribution of the relative momentum k = 
(ki -A:2)/2 of the particle pair with total spin S, Ms and isospin 
r, Mt is defined as 
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We also define the two-body densities p^g\{r) and rffj{k) that 
are obtained by summing the spin- and isospin-indices Ms 
and Mt so that they do not depend on the orientation of r and 
k 
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The distributions defined in Eqs. (7) and (8), when cou- 
pled properly in space-spin space, are called internucleon cor- 
relation functions in Ref. [25]. The internucleon correlation 
functions contain all the information needed to calculate the 
energy of the state for a two-body Hamiltonian. 



C. Realistic nuclear forces 
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where p^ ^ y,^ (^i , ^2 ) is the probability density that a nucleon 
pair with one nucleon at position ri and the other one at r2 is 
found in the spin S, Ms and isospin T, Mt channel, ri and r2 
are measured from the total center-of-mass position and P.!^^ 

and p™^ project on spin and isospin of the pair, respectively. 
The label T here indicates the two-nucleon isospin. Note that 
it is also used to denote the total isospin of the system in the 
previous subsection. 

To keep the graphical presentation transparent we discuss 
the short-range correlations as a function of the relative posi- 
tion r = ri - r2 of the two nucleons only and integrate over 



The Argonne v8' (AYS') potential [1] is depicted in Fig. 1 
as a function of r for the four spin-isospin combinations of 
a nucleon pair. In the left most graph for 5" = 1 , T = we 
assume the nucleons to be at rest so that the spin-orbit in- 
teraction does not contribute. The tensor interaction causes 
a quadrupole type dependence as a function of the angle be- 
tween the total spin direction (which we aligned along the z 
axis) and the direction of the distance vector r. The main at- 
traction is obtained when the spins of the nucleons are aligned 
with the distance vector r while almost no attraction exists 
in the x direction where the spins are orthogonal to r. For 
S = 1 , r = 1 we added the spin-orbit interaction for L^ = I 
because due to the Pauli principle the nucleon pair has to be 
in an orbital state with odd parity. For T = 1 the tensor inter- 
action leads to a situation which is opposite to the T = case. 
Here the attraction occurs along the x axis where the spins are 
orthogonal to r. 

In the 5" =0 channels there is no tensor interaction and 
no spin-orbit interaction thus the interaction depends only on 
the distance \r\. Common to all channels is the strong central 
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FIG. 1. (Color online). Argonne v8' potential in the different spin-isospin channels as a function of the distance vector r = (x, j = 0,z). In the 
5 = 1 channels the total spin is aligned with the z axis. Units are in MeV. 



repulsion for \r\ < 0.6 fm. For S = 0,T = 1 there is strong 
attraction, around \r\ = 1 fm, however, not strong enough to 
make the di-neutron bound. The S = 0, T = potential is 
repulsive at all distances. 

This rather complex nature of the nuclear interaction in- 
duces corresponding intricate correlations in the A-body 
eigenstate of the Hamiltonian which we discuss in Sec. III. 



III. RESULTS 

In the following we investigate the ground states of ^H with 
r = 1+, ^H and ^He with r = |^, and ^He with r = 0+, 
labeled by d, t, h, and a, respectively, as well as the excited 0^ 
state of "^He at 20.21 MeV, which is a resonance close to the 
threshold for ^H -\- p, labeled by a*. In this paper we treat the 
state a* as a quasi-bound state with a long tail [31], though it 
has a proton width of 0.5 MeV. 



A. One-body densities 

The one-body point densities of the five states are depicted 
in Fig. 2. In all cases the position ri of the nucleon is counted 
from the total center-of-mass position of the many-body sys- 
tem. For the deuteron this means that ri = jr is half the 
relative distance between neutron and proton. The densities 
are averaged over the directions of the total spin and hence 
depend only on ri = |ri|. Likewise the momentum ki of a 
nucleon is the one in the total center of momentum frame and 
averaging over total spin directions is implied. 

The a particle shows the largest central density, the ^H and 
^He densities are somewhat smaller and diff'er only slightly 
due to the Coulomb interaction. The density of the excited 
0^ state in "^He is much lower because this state, which is a 
narrow resonance in the scattering of ^H and proton, is es- 
sentially a configuration in which a proton and a triton orbit 
around each other in an / = state [3 1,32]. Due to the recoil 
the quantal zero point motion in the relative coordinate smears 
out the intrinsic density of the triton. 
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FIG. 2. (Color online). One-body point densities of the different 
states in coordinate space (top) and one-body density in momentum 
space divided by mass number A (bottom). Ground states of ^H, ^H, 
^He, ^He are denoted by d, t, h, and a, respectively. The excited state 
of "^He is labeled with a*. 



We include also the deuteron despite the fact that its one- 
body density is actually the two-body density at half the dis- 
tance, and only the 5" = 1,7 = component of the four pos- 
sibilities to couple spins and isospins of two nucleons con- 
tributes. The comparison with the three- and four-body sys- 
tems nicely demonstrates that in coordinate space the eff'ects 
of the short-range repulsion, which are clearly visible in the 
deuteron, can not be seen in the A-body case because the one- 
body density integrates over the positions of the other A - I 



particles. 

However the one-body momentum distribution (lower part 
of Fig. 2) shows beyond ^i ^ 1.5 fm"^ the presence of short- 
range correlations by a far-out-reaching tail. The occurrence 
of high momenta is at variance with a Hartree-Fock like 
mean-field picture where beyond the Fermi momentum kp ^ 
1.4 fm"^ the momentum distribution drops steeply [18, 33]. 
For the deuteron the two-body density is again identical to the 
one-body density (in momentum space ki - -k2 = k). One 
notices that the high momentum tails have a very similar form 
in all cases including the deuteron. This similarity suggests 
already a universal behavior of the short-range correlations 
independent of the spatial density of the A-body system. 



B. Two-body densities 

The A-body density, which contains the information about 
all correlations, is a function of A position or momentum vec- 
tors and 4A spin-isospin possibilities and hence can not be 
visualized easily. Therefore we integrate and sum over A - 2 
single-particle degrees of freedom and are left with the two- 
body density. This represents the sum over all particle pairs in 
the many-body state. In addition we integrate over the center- 
of-mass position of the pair and obtain the two-body densi- 
ties defined in Eqs. (7) and (8) for the four spin-isospin chan- 
nels which are possible for a nucleon pair. The complex na- 
ture of the nucleon-nucleon interaction discussed in Sec. II C 
induces short-range repulsive and tensor correlations in the 
many-body state, which can be seen best in the two-body den- 
sity. 

In Fig. 3 the spatial two-body densities P^iIqq(t) of the four 
diff'erent states are displayed. The first striking observation is 
that at short distances they look very similar independently of 
the many-body state. That means that the correlations felt by a 
particle pair in the 5" = 1,7 = channel are at short distances 
the same independently of the remaining particles in the sys- 
tem. The second not unexpected observation is that these den- 
sities reflect in an almost one-to-one fashion the potential in 
the 5" = 1, r = channel, see Fig. 1. There exists a one-to- 
one correspondence between the nuclear Hamiltonian and the 
two-body densities. The expectation value of the Hamiltonian 
can be calculated with the two-body density as discussed in 
Ref. [25]. In regions where the potential is most attractive, 
r ^ (0, 0, ±1 fm), the densities are large and where the inter- 
action is repulsive or close to zero the probability of finding 
the particle pair is small. At small distances below 0.5 fm 
the AV8' potential is so strongly repulsive that the pair den- 
sities in all many-body states are pushed down toward zero. 
One should bear in mind that in a simple shell model many- 
body state these correlations can not be represented. The shell 
model two-body densities have actually their maximum at rel- 
ative distance r = 0. 

For 5" = 1 , M^ = 1 , r = 1 the tensor interaction is most at- 
tractive in a torus around the z axis (see Fig. 1) and hence the 
two-body density has its maximum in the x-y plane as can be 
seen in Fig. 4. In this channel we see again a one-to-one cor- 
respondence to the potential. For small distances up to about 



TABLE L Number of pairs in different states of light nuclei calcu- 
lated with the AV8' potential. Calculated binding energy E^ in MeV, 
matter point radius ^j{r^) in fm. 
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1.5 fm the shape of the distribution is again almost identi- 
cal for all three many-body states. It should be noted that 
this channel, which does not exist for the deuteron, is weakly 
populated in the three- and four-body systems. Depending on 
the nucleus about 5-7% of the pairs are in this channel, see 
Table I. In the shell model representation this channel corre- 
sponds to at least one particle-hole excitation to the p shell 
such that the relative motion of the pair has negative parity. 
When occupying only the s shell this channel is forbidden. 
We will discuss in Sec. HIE how the two-body density in this 
channel is related to three-body correlations. 

The second strong channel has S = 0, T = 1 and is shown 
in Fig. 5. As there is no tensor and spin-orbit interaction for 
S = the distributions are spherical. Again they are very 
similar for all states and also exhibit a hole at short distances 
where the AV8' potential is very repulsive and a maximum at 
distances around 1 fm where it is most attractive, see Fig. 1. 

The S = 0, r = channel is not displayed because its 
contribution listed in Table I is tiny, only about 0.1%. The 
potential in this channel is purely repulsive as can be seen 
in Fig. 1. Nevertheless this small contribution is surprising 
when compared with the 5" = l,r = 1 channel where the 
potential, while not purely repulsive, provides only very weak 
attraction. We will discuss this point later in relation to three- 
body correlations. 

It is also interesting to note that the number of pairs in the 
O2 state of the ^He nucleus are almost identical to the summed 
number of pairs from the triton and ^He — reflecting the clus- 
ter nature of this state. 

In Figs. 6 and 7 we show the two-body densities in momen- 
tum space in the 5 = 1, T = and 5" = 0, T = 1 channels. As 
expected we find more pairs at small relative momenta for the 
spatially extended deuteron and the excited 0^ state in ^He. In 
all states we observe extended high momentum components 
above the Fermi momentum of about 1.4 fm"^ Comparing 
both even channels we find that the two-body densities are 
very similar at low momenta up to about 0.5 fm"^ but the 
high-momentum components in the 5" = 1 , T = channel 
are larger by a factor of 2-2.5. The diff'erences in the number 
of pairs (see Table I) in the even channels originates mainly 
from contributions at higher momenta between 0.5 and about 
2.5 fm"^. The larger number of high momentum pairs in the 
S = 1 , r = channel can be traced back to tensor correlations 
as we will discuss in Sec. IIIF. 

These diff'erences in the high-momentum contributions be- 
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FIG. 3. (Color online). From left to right: Two-body densities in coordinate space for a pair of nucleons with 5 = 1, M5 = 1 and T = in the 
ground states of ^H, ^H and "^He and the 20.21 MeV excited state of "^He denoted by d, t, a, and a*, respectively. The densities have rotational 
symmetry around the z axis and range from black = to bright (yellow) - maximum. Maxima assume values of 0.008 fm~^ for d, 0.015 fm~^ 
for t, 0.035 fm-3 for a, and 0.015 fm'^ for a*. 




FIG. 4. (Color online). From left to right: Two-body densities in coordinate space for a pair of nucleons with *S = 1, Ms - 1 and T = 1 of ^H 
and "^He in the ground states and the 20.21 MeV excited state of ^He denoted by t, a, and a* respectively. The densities are axially symmetric 
around the z axis and Mt is summed over. Otherwise same as Fig. 3. Maximum densities are 0.47 x 10~^ fm~^ for t, 2.2 x 10~^ fm~^ for a, and 
0.51 X 10"^ fm-3forcK*. 




FIG. 5. (Color online). The same as Fig. 4 but for a pair of nucleons with 5=0 and T - I. Maximum densities are 0.020 fm ^ for t, 
0.054 fm-3 for a, and 0.019 fm'^ for a*. 




TABLE 11. Normalization factors C^j. 



FIG. 6. (Color online). Two-body densities as a function of relative 
momentum k for the S = l,T = channel. Ground state densities of 
^H, ^H, ^He, ^He are denoted by d, t, h, a, respectively. The excited 
state of ^He is labeled with a*. 




FIG. 7. (Color online). The same as Fig. 6 but for the 5" = 0, T = 1 
channel as a function of k. 



tween the 5" = 1,7 = and S = 0, T = 1 channels can also be 
interpreted in terms of proton-proton or neutron-neutron (only 
r = 1) versus proton-neutron pairs (both T = and 7 = 1). 
Such effects have been found also in theoretical studies for 
heavier nuclei [19, 22, 23, 34] and in experiment comparing 
for example the (e, e'pp) with (e, e'pn) cross sections where a 
dominance of proton-neutron pairs was observed [11]. 



C. Universality at small distances 

As already seen in Figs. 3-5 the two-body densities of 
the different states look very similar especially at small dis- 
tances. To further investigate this universality of the short- 
range correlations we display in Fig. 8 cuts of the two-body 
density P^ilo^ir) along the z- and the x-direction. As the ab- 
solute values of the densities are quite different in the five 
states (see Fig. 3) we normalize the two-body densities at 
r = 1 fm, where the densities approximately reach their max- 
imum value. The normalization factors 
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FIG. 8. (Color online). Cuts of the normalized densities pf/oo^'') ^^^ 
r = (x,0,0)andr = (0,0, z). 



are given in Table II. This choice for the normalization radius 
is not crucial as the ratios of the normalization coefficients be- 
tween different states for a given channel are essentially con- 
stant (within 2%) when calculated between and 1.0 fm. 

It is astonishing to see in Fig. 8 that for small distances the 
scaled densities practically coincide along both the z and the x 
axes. This means that not only the central correlations but also 
the angular-dependence of the tensor correlations are almost 
identical at short distances. The short-range central and ten- 
sor correlations exhibit universal behavior at short distances 
below about 1 fm. 

In the 5" = 0, r = 1 channel the same universal behavior 
can be observed as shown in Fig. 9. The two-body densities 
normalized at r = 1 f m for the different systems agree per- 
fectly to distances up to about 1 fm. 

Whereas the behavior of the two-body densities at short dis- 




FIG. 9. (Color online). Two-body densities Po^l(r) normalized to 
1 fm~^ at r-l fm for different states (c.f. Fig. 5). 




FIG. 10. (Color online). Normalized two-body densities as a func- 
tion of relative momentum k for the S = l,T = channel. Ground 
state densities of ^H, ^H, ^He, ^He are denoted by d, t, h, a, respec- 
tively. The excited state of ^He is labeled with a*. 




FIG. 11. (Color online). The same as Fig. 10 but for the 5" = 0, T = 1 
channel as a function of k. 



tances is universal the behavior at large distances is specific 
for the particular many-body state. Its form is discussed in 
Ref. [25]. 



D. Universality at high momenta 

In Fig. 10 and Fig. 11 we show the two-body densities in 
momentum space in the 5" = 1,T = and S = 0,T = 1 
channels scaled with the same normalization factors as given 
in Table II that were determined for the two-body densities in 
coordinate space. Whereas the scaled densities diff'er strongly 
at low momenta, we find almost perfect agreement at high mo- 
menta larger than about 3 fm"^ The universality of the short- 
range correlations in coordinate space is therefore reflected in 
a universality of the high-momentum components in momen- 
tum space. The fact that the two-body densities difl'er in the 
intermediate momentum range from the Fermi momentum of 
about 1 .4 fm"^ to about 3 fm"^ should be related to difl'erences 
in the long-range correlations for the difl'erent nuclei. 



E. Three-body correlations 

Looking at the number of pairs in the difl'erent spin-isospin 
channels (see Table I) an interesting observation can be made. 
Let us concentrate on "^He. In a simple shell model picture 
where all nucleons occupy s orbits we should find three pairs 
each in the 5" = 1 , r = and S = 0,T = 1 channels and 
zero pairs in the odd channels. The nuclear potential is much 
more attractive in the even channels than in the odd chan- 
nels, furthermore the kinetic energy is much higher in the odd 
channels due to the non- vanishing angular momentum. It is 
therefore surprising that we find in the exact wave function a 
remarkable depopulation of the 5" = 0, T = 1 even channel 
(2.572 pairs) obviously in favor of the 5" = 1, T = 1 odd chan- 
nel (0.428 pairs). As remarkable is the fact that the number 
of pairs in the 5" = 1, T = channel is essentially unchanged 
(2.992 pairs) compared to the simple shell model picture. This 
eff'ect can not be understood in terms of two-body correla- 
tions, as the parity of the relative motion of a nucleon pair can 
not be changed by the two-body interaction. As already dis- 
cussed by Forest et al [21] this eff'ect should be attributed to 
three-body correlations induced by the strong tensor force in 
the 5" = 1, r = channel. As total isospin T is a conserved 
quantity in light nuclei the total number of pairs in the T = 
and r = 1 channels has to be conserved. The tensor force in 
the 5" = 1 , r = channel provides the dominant contribution 
to the nuclear binding. It has its origin in the pion exchange 
and is long ranged. Nucleon pairs in the 5" = 1,7 = channel 
will therefore be correlated even at large distances and these 
tensor correlations will aff'ect other nucleon pairs. It is ener- 
getically favorable to break a pair in the 5" = 0, T = 1 channel 
by flipping the spin of a nucleon if this allows the tensor force 
to gain energy in another pair involving a third nucleon. An 
illustration of this mechanism is shown in Fig. 12 where en- 
ergy is gained by tensor correlations for a pair of nucleons 
in the S = 1,T = channel. In the uncorrelated case the 
nucleon pair is assumed to be in a relative S -wave. In the cor- 
related many-body state the pair will be partially found in a 
relative Z)-wave to allow for additional binding by the tensor 
force. This Z)-wave admixture will also change the spin ori- 
entation of the nucleons, so that another pair, originally in the 
S = 0, r = 1 channel, is now found in the S = 1 , T = 1 
channel. 

To illustrate the efl'ects of these three-body correlation on 
the two-body densities in the T = 1 channel we show in 
Fig. 13 the two-body momentum distributions of the S = 
0, r = 1 and the 5" = 1, r = 1 channels for "^He. At smaU 
relative momenta the density in the odd channel vanishes be- 
cause of the P-wave nature. For momenta between 1.1 and 
2.1 fm"^ the two-body density in the 5" = 1, T = 1 is actually 
larger than in the 5" = 0, T = 1 channel. At very high rel- 
ative momenta the contribution of the odd channel can again 
be neglected. The three-body correlations therefore influence 
the two-body density very difl'erently in difl'erent momentum 
regimes. For low relative momenta below about 0.5 fm"^ the 
efl'ect is very small and the two-body densities in the two 
even channels are very similar. In an intermediate momen- 
tum range between 0.5 and 2.5 fm"^ we observe a notice- 
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FIG. 12. (Color online). Illustration of three-body correlations in- 
duced by tensor correlations. In the uncorrelated wave function (left) 
the two nucleons 1 and 2 are in an *S = 1 , Ms - pair with L = 0. 
The tensor force leads to an admixture of an L = 2 component and 
an alignment of the spins of nucleons 1 and 2 flipping the spin of 
nucleon 2 (right). This affects the interaction between nucleon 2 and 
nucleon 3. In the uncorrelated wave function the protons 2 and 3 
form an*S=0, r = l,L = pair. After the spin-flip of nucleon 2 
this becomes miS = l,r = 1,L= 1 pair. 



body correlations. Smaller three-body correlations means that 
less nucleon pairs are moved from the 5" = 0,T = 1 to the 
S = 1, r = 1 channel. As in the odd channel the potential 
is less attractive and the kinetic energy is much larger, the 
three-body correlations provide a repulsive contribution to the 
energy. 

It has already been realized that a term in the effective inter- 
action called antisymmetric spin-orbit (ALS)-force that con- 
nects S = with S = I states and changes the relative angu- 
lar momentum by AL=1, like (h - h) • (o"i - 0-2), is able to 
improve spectra and transition rates in ^J- shell model calcu- 
lations [40] . But as such a term is not conserving translational 
and Galilei invariance it is not allowed in the free nucleon- 
nucleon interaction and can only be obtained by integrating 
many-body forces over additional particle degrees of freedom. 

We want to stress the point that in our discussion no gen- 
uine three-body forces are considered. The three-body cor- 
relations are induced by the two-body tensor force. When 
genuine three-body forces are included we of course expect 
additional or modified three-body correlations. 




FIG. 13. (Color online). Two-body densities in momentum space for 
"^He in the iS = 0, T = 1 and S - \,T - 1 channels and the sum of 
both densities. 



R Comparison with unitary correlation operator method 

The universality of short-range correlations is not only in- 
teresting in itself but also confirms the basic assumptions that 
underlie methods to derive effective low-momentum interac- 
tions such as UCOM, Viow-k and SRG. We will discuss here 
the UCOM approach as it provides the most direct connection 
to the short-range correlations in the nucleus. 

The basic idea of the UCOM approach is to imprint the 
short-range central and tensor correlations into the nuclear 
many-body wave functions explicitly by means of a unitary 
correlation operator C. Starting from an uncorrelated trial 
state |0) the correlated state 



|^) = C|0) 



(12) 



able depletion of the 5" = 0, T = 1 channel in favor of the 
S = \,T = \ channel. This contributes to the fact that the 
two-body densities in the 5" = 1 , T = channel are much 
larger than in the 5" = 0, T = 1 channel in this momentum 
region. 

As already mentioned this effect can not be understood in 
terms of two-body correlations. It also explains why effective 
interactions that are obtained by unitary transformations in 
two-body approximation, like Viow-k [35], the similarity renor- 
malization group (SRG) [36] or the unitary correlation opera- 
tor method (UCOM) [37-39], provide more binding than the 
bare interaction when used in exact calculations. In two-body 
approximation the interaction is transformed independently in 
all spin-isospin channels. It is therefore possible to obtain the 
full contribution of the tensor force in the S =1,7 = chan- 
nel without having to pay the price of the three-body corre- 
lations. With increasing range of the tensor correlations (in 
the UCOM approach) or a lower cut-off (in the Viow-k or SRG 
approaches) the effective interaction will induce smaller three- 



then features the short-range central and tensor correlations. 
Long-range correlations still have to be incorporated explic- 
itly in the trial state |0). 

To explain the action of the correlation operators we discuss 
first how the relative motion of two nucleons is affected by the 
correlation operators. For that we use basis states 



\(P(LS)JM;TMt}. 



(13) 



where the relative orbital angular momentum L is coupled 
with the spin S of the two nucleons to total angular momen- 
tum /, M. The isospin is coupled to T, Mr. The radial part of 
the relative wave function is given by (p(r). 

In the 5" =0 channels only the central correlation operator 
acts and the correlated relative wave function is given, using 
the correlation function R-(r), as 



^i'' 



(r) = {riLS)JT\ C, \</){LS)JT) 



(14) 
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whereas in the 5" = 1 channels both central and tensor corre- 
lation operators act and we obtain the correlated radial wave 
functions 



fuuir) = {r(L'S)JT\ C^Cr \(f>(LS)JT} = 

;U =L = J±l 

\L = J±l,L=J+l 

(15) 



0.08 



cosOj(r)^^lRUF)cp(R.(r)) 
±smOj(r)^ ^/RUr) (P(R.(r)) 



with the tensor correlation function 



6J(r) = 3^fJ(JTT)§(r). 



(16) 



The functions R-(r) and ?^(r) also carry implicitly the ap- 
propriate quantum numbers which are omitted here. 

To calculate the two-body density in momentum space we 
will need the relative wave function in momentum space as 
obtained by Fourier transformation 



{q(VS)JT\CnCr\(p(LS)JT) = 

drr'ju{qr){r{L'S)JT\ CnCr \^{LS)JT) . 



#T 



(17) 



To illustrate the action of the correlation operators we re- 
strict the discussion here to the most simple trial state for "^He 
where all nucleons occupy the s orbit in a harmonic oscillator 



m = 1(05)4) _ 



(18) 



The harmonic oscillator width parameter a = 1.98 fm is ad- 
justed to reproduce the radius of the "^He nucleus as obtained 
in the exact calculation with the correlated Gaussian approach. 
We can then express the uncorrected two-body density op- 
erator for this state as 



/^uncorr 



M 



with the relative L ■ 
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Including short-range central and tensor correlations with 
the UCOM correlation operators the two-body density opera- 
tor of the correlated state is given in two-body approximation 
as 



a(2) _ /^ /^ ^K^) /^T/^ 



(2) r^M 

uncorr ^r^Q 



Y, l<Ao'"(00)0; lMr)<.Ar'(00)0; IMt\ + 

Mj 
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(21) 
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FIG. 14. (Color online). Comparison of the coordinate space two- 
body density in the 5" = 0, T = 1 channel in "^He between the UCOM 
and the exact many-body calculation using correlated Gaussians de- 
noted by CG See the text for the different UCOM and uncorrected 
results. 
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FIG. 15. (Color online). The same as Fig. 14 but for the two-body 
density in the »S = 1,7 = channel as a function of x and z. 



where L = 2 components appear in the S = 1,7 = channel 
due to the tensor correlation operator. 

In the following we use correlation functions derived from 
an SRG evolved AV8' Hamiltonian [39, 41]. The used flow 
parameters a = 0.04 fm"^ and a = 0.20 fm"^ correspond to cut- 
ofl' parameters of /I ^ 2.2 fm"^ (soft) and A ^ 1.5 fm"^ (very 
soft). The labels UCOM04 and UCOM20 will be used in the 
following to identify these two sets of correlation functions. In 
exact calculations using the no-core shell model [39] the cor- 
responding UCOM interactions provide binding energies that 
are close to the experimental binding energies for -^H and "^He. 
Using the simple trial state in Eq. (18) we obtain "^He binding 
energies of -18.50 MeV and -25.10 MeV with UCOM04 and 
UCOM20, respectively. 

In Figs. 14 and 15 we compare for the two even channels 
the two-body densities in coordinate space given by Eq. (21) 
with the two-body densities calculated from the exact solution 
for ^He. The two-body densities obtained within the UCOM 
approach agree very well with the exact two-body densities at 
small distances. Compared with the uncorrected wave func- 
tion the two-body density is strongly suppressed at short dis- 
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FIG. 16. (Color online). UCOM two-body density in momentum 
space for the 5" = 0,T = 1 channel in "^He compared with exact 
many-body calculation denoted by CG. 
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FIG. 17. (Color online). The same as Fig. 16 but for the *S = 1,7 = 
channel. Contributions from the L = and L = 2 components are 
shown for the UCOM densities. 



tances. At r around 1 fm the UCOM two-body densities de- 
pend on the specific choice of the correlation functions and 
deviate by 0-20% from the exact results. In the 5 = 0, T = 1 
channel the UCOM two-body densities are always larger than 
the exact results. The main reason for this discrepancy is the 
two-body approximation that we used to obtain the correlated 
two-body density. The UCOM results do not include the ef- 
fects of three-body correlations and the number of pairs in the 
presented UCOM result is exactly three in both even channels. 

In the S = 1,T = channel (Fig. 15) we find an al- 
most perfect agreement of the exact two-body densities with 
the UCOM20 result which uses long-range correlation func- 
tions. This holds not only for the radial dependence due to 
short-range repulsion but also for the angular dependence of 
the two-body density due to the tensor correlations. For the 
shorter-ranged correlation functions UCOM04 the agreement 
is not so good. This is caused by the diff'erent ranges of the 
tensor correlation functions. It appears that the long-range 
correlation functions in UCOM20 are able to describe most of 
the tensor correlations found in ^He whereas with the short- 
ranged tensor correlation functions in UCOM04 a significant 
part of the medium to long-range tensor correlations is miss- 
ing. 

The two-body densities in momentum space obtained with 
the UCOM densities given in Eq. (21) are compared with the 
exact results in Fig. 16 and Fig. 17. The eff'ect of the short- 
range correlations are reflected in the high-momentum com- 
ponents. For relative momenta larger than about 3 fm"^ we 
find good agreement with the exact result for both UCOM 
transformations. In the 5" = 0, T = 1 channel the UCOM 
densities are much too small in the intermediate momentum 
region from 1.4 to 3 fm"^ This reflects the oversimplified 
Gaussian trial wave function and missing many-body correla- 
tions. 

For the two-body densities in the 5" = l,T = chan- 
nel (Fig. 17) we have decomposed the UCOM results in the 
L = and the L = 2 components. The L = component 
looks very similar to the two-body densities obtained in the 
S = 0,T = I channel and does not contribute in the interme- 
diate momentum-region. The L = 2 component, introduced 



by the tensor correlation operator, on the other hand actually 
dominates the medium- and high-momentum part of the two- 
body density. There is a strong dependence on the range of 
the correlation functions but even in the UCOM20 case the 
exact two-body densities are still significantly larger in the 
intermediate momentum region. Again, contributions due to 
long-range correlations are missing. 

To include these missing contributions consistently the 
UCOM two-body densities should be calculated not from the 
simple trial state in Eq. (18) but from an exact solution |0) of 
the many-body problem 



HucoMm = Em 



(22) 



using the UCOM eff'ective Hamiltonian in two-body approxi- 
mation 



^ucoM = C HC . 



(23) 



Such calculations are in preparation using the no-core shell 
model. Nevertheless we can not expect perfect agreement 
even in this case due to the two-body approximation. Using an 
eff'ective interaction like UCOM the three-body correlations as 
discussed in Sec. HIE will not be fully included. The devia- 
tions between the exact two-body densities and the two-body 
densities obtained using eff'ective interactions in two-body ap- 
proximations will depend on the range of the correlation op- 
erators (in the UCOM approach) or on the value of the cut-off* 
(in SRG and Vw-yi). 

A detailed discussion of the operator evolution in the SRG 
for the deuteron is provided by Anderson et al. [42]. The 
authors study the evolution of high-momentum operators like 
the momentum distribution within the SRG and investigate to 
what extend a decoupling between low- and high-momentum 
components occurs. 



IV. CONCLUSIONS 

In this paper we have studied the two-body densities in co- 
ordinate and momentum space for the deuteron, ^H, ^He, ^He 
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and the first excited 0"^ state in "^He. Fully converged solu- 
tions for these light nuclei could be achieved using the corre- 
lated Gaussian basis approach for the Argonne v8' interaction. 
The short-range repulsion and the tensor force induce strong 
short-range correlations in the many-body wave functions, re- 
flected in the two-body densities. If the two-body densities 
in coordinate space are normalized at short distances, we find 
in the diff'erent spin-isospin channels a universal behavior up 
to about 1 fm in all nuclei. Using the same normalization we 
observe a corresponding universal behavior of the two-body 
densities in momentum space at relative momenta larger than 
about 3 fm"^. Although we only have two-body forces we 
could identify three-body correlations due to the long-range 
tensor correlations in the S = 1,7 = channel. They man- 
ifest themselves in the two-body densities by a reshuffling of 
pairs from the 5" =0,7= 1 channel into the 5 = 1 , T = 1 
channel. 

The universal behavior of the short-range correlations ex- 
plains the success of approaches such as Viow-k [35], or SRG 
[36] and UCOM [39] that use unitary transformations to de- 
rive an eff'ective low-momentum interaction. The idea of the 
unitary transformation is to decouple the short-range from the 
long-range or the high-momentum from the low-momentum 
physics. Using such transformed low-momentum interactions 
the wave functions no longer show the strong short-range cor- 
relations induced by the original interaction. To recover short- 
range correlations the two-body densities have to be trans- 
formed using the same unitary transformation. We compared 
in coordinate and momentum space the exact two-body den- 
sities of "^He to those obtained from a simple OHQ. trial wave 
function and the UCOM transformation for the Argonne v8' 
interaction. In the S = 1,7 = channel we find a very good 
agreement for the short-range and the high-momentum behav- 
ior of the two-body densities. Diff'erences show up mostly in 
the intermediate momentum range from 1.5 to 3 fm"^ In this 
region long-range correlations, missing in the very simple trial 
wave function, become important. In the 5" = 0, T = 1 chan- 
nel the agreement is spoiled by the missing three-body corre- 
lations in the UCOM approach. To recover these diff'erences 
the unitary transformation would have to be performed not in 
a two-body approximation as done here but on the three-body 



level. 

In a more elaborate approach short-range correlations in 
heavier nuclei will be studied by solving the many-body prob- 
lem with a soft unitarily transformed interaction for example 
with the no-core shell model [39] and then calculating the 
unitarily transformed two-body densities. This will allow to 
include both, long-range correlations by the many-body ap- 
proach and short-range correlations by the unitary transfor- 
mation. In the two-body approximation the role of three-body 
correlations could be investigated by varying the cut-off' of the 
transformation. An explicit treatment of three-body correla- 
tions is possible in principle but would become very involved. 

We studied the two-body densities in this paper only as a 
function of the distance or the relative momentum of the nu- 
cleons, but it would also be interesting to investigate the de- 
pendence on the center of mass momentum of the nucleon 
pairs. Wiringa et ah found that the short-range correlation 
eff'ects are most pronounced at vanishing center-of-mass mo- 
mentum for the pairs [23]. For larger center-of-mass momen- 
tum the short-range correlations are smeared out as there is a 
higher probability to find one nucleon inside the Fermi sphere 
even at high relative momentum. It might also be interesting 
to study the two-body densities as a function of the center- 
of-mass position. In particular it might be possible to study 
short-range correlations of neutrons in the surface of neutron- 
rich exotic nuclei. Because of the universality of the short- 
range correlations information from this low-density regime 
should also be important for the saturation properties of neu- 
tron matter at higher densities. Of course three-body forces 
will become more and more important with increasing den- 
sity. 
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